home *** CD-ROM | disk | FTP | other *** search
/ Magnum One / Magnum One (Mid-American Digital) (Disc Manufacturing).iso / d18 / nrpas13.arc / QGAUS.DEM < prev    next >
Text File  |  1991-05-01  |  537b  |  31 lines

  1. PROGRAM d4r7(input,output);
  2. (* driver for routine QGAUS *)
  3. CONST
  4.    x1=0.0;
  5.    x2=5.0;
  6.    nval=10;
  7. VAR
  8.    dx,ss,x : real;
  9.    i : integer;
  10.  
  11. FUNCTION func(x: real): real;
  12. BEGIN
  13.    func := x*exp(-x)
  14. END;
  15.  
  16. (*$I MODFILE.PAS *)
  17. (*$I QGAUS.PAS *)
  18.  
  19. BEGIN
  20.    dx := (x2-x1)/nval;
  21.    writeln;
  22.    writeln('0.0 to','qgaus':10,'expected':13);
  23.    writeln;
  24.    FOR i := 1 to nval DO BEGIN
  25.       x := x1+i*dx;
  26.       qgaus(x1,x,ss);
  27.       writeln(x:5:2,ss:12:6,(-(1.0+x)*exp(-x)
  28.          +(1.0+x1)*exp(-x1)):12:6)
  29.    END
  30. END.
  31.